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O ! Abstract 

^H ' We analysed the identified hadron multiplicities predictions of the generalised thermodynamical 

■^ ■ model of the multiparticle production processes with non-extensive statistic. The multiplicities 



measured recently at LHC experiments are found to be consistent with this approach and ther- 



,^ \ modynamical parameter values found already analysing transverse momentum distributions. The 

Ph. 
jiv ' information about the mechanism of the strangeness suppression has been derived to some extend 

(D ■ 

r"! ■ from multistrange hadron abundances. 
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INTRODUCTION 

The thermodynamical picture of particle creation in hadronic collisions was the first and 
quite successful attempt to describe it. The elaborated and complete theory of this kind 
was presented in serial of papers by Hagedorn [see [1] and references therein]. The idea of 
the fireball with the proposition that "all fireballs are equal" give very strong predictions 
concerning produced particle spectra. 

One of the predictions was that the temperature of the "hadronic soup" (precisely defined) 
could not exceed a universal constant Tq of order of 160 MeV. This value comes from exam- 
ination of elementary particle mass spectrum not as a result of the procedure of parameter 
adjusting using multiparticle production (e.g., transverse momenta) data. The Hagedorn 
theory were abundant (for some time) when more sophisticated, jet or QCD based ideas 
appeared [2]. One of the reasons was the failure of the high-transverse momenta descrip- 
tion. The temperature of the fireball is defined as the parameter in the classical Boltzmann 
exponential term of the probability weights for phase space average occupation numbers. 
This defines the (asymptotic) form of the distribution of transverse momentum of particles 
created from decaying fireballs. It was found that at high and very high interaction energies 
the predicted exponential fall do not agree with the observed high p^ behaviour. 

Successes of QCD based description of the hard processes gave deep insight into the 
nature of physics involved, and common believe that this is just the right theory of strong 
interactions, makes the thermodynamical approach only one, very approximate, simple and 
naive, tool of limited applicability and thus limited significance. But on the other hand, 
the simplicity of the theory and notorious constant lack of the effective QCD theory of soft 
hadronization processes give a hope that the fireball idea can be enriched, modified and can 
become important again. 

The identified hadron ratios have been maesured with all LHC detectors and results were 
compared with high-energy event generators available in the market. The comparison was 
not very satisfactory. 

The ATLAS results on strange particle production are compared with various Monte 
Carlo Simula- tion models [3] and again some substantial disagreements are found in distri- 
butions of transverse momentum. 

The LHCb prompt hadron measurement ratios [4] were compared with several Monte 



Carlo generator results and it was concluded that none was able to describe well all observ- 
ables, e.g., significant underestimation of {K'^ + K^) / [ti^ + 7r~) has been observed with the 
Perugia and Perugia NOCR generators. 

The ALICE data on K* ans ip production [5] have been compared to a number of PYTHIA 
[6] tunes and the PH0JET[7] event generator. None of them gives a fully satisfactory 
description of the data. 

ALICE Collaboration measured also neutral pion and -q mesons [8] and concluded that 
calculations cannot provide a consistent description at 7 TeV, as well multi-strange baryons 
Ss and Vis [9], and found that particle yields, (j5_|_), and the spectra in the intermediate p^ 
range are not well described by the PYTHIA Perugia 2011 tune Monte Carlo event generator 
specially for Vt. 

PARTITION FUNCTION 

The Hagedorn idea was used again to describe the identified particle multiplicities in 
hadronization both, in e^e~ annihilation and hadronic collisions. The grand canonical for- 
malism of Hagedorn was replaced in the serial of papers by Becattini and co-workers [10] by 
the canonical one, much relevant for studies of small systems like primary created fireballs 
for which the requirement of exact conservation of some quantum numbers seems important. 
In general, thermodynamics of the system is determined by the partition function which can 
be written as 

Q i,j 

where x is the Boltzmann exponential factor x = exp {—Ej^/T), j and k enumerate particle 
types and momentum cells, Q^ is the initial fireball quantum number vector and Q is the 
respective vector of the particular state, and Uj^ is the occupation number. For bosons 
one has Ujk = 0, 1, 2, 3 . . . and ^ x/^ = 1/(1 — Xjk), while for fermions Ujk = 0, 1 and 

J k J 

J2,y- ^]k = 1 + ^ik- Thus introducing Fourier transform of b (and reducing vector Q to 



3-dimensional: charge, baryon number and strangeness) Eq.(l) becomes 
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where Qj is the quantum number vector of the particle j. 

The next step is to turn the summation over phase space cells to an integration over 
momentum space 

Y.^-.j^J<^'p (3) 

where Wj is the weight factor associated with the particle j. The first guess is that it should be 
equal to (2Jj + l) and counts spin states. However, this does not to be so simple (see, e.g., [6]) 
and other solutions introducing factors responsible for some wave-function normalisation, 
which should disfavour heavier states were found to be preferable by measurements. We will 
discuss this point later on. 

With the Eq.(2) we are ready for detailed numerical calculations. 

Average multiplicities - Boltzmann statistic 

With the known partition function Z the average characteristics of the system can be 
obtained in an usual way. The average occupation number Uj representing the mean multi- 
plicity of the particles of the kind j is (introducing the fictions fugacity A) given by 



(nj) = |^log(Z(Q°, A)) 



(4) 

A=l 



thus 
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where the upper sign is for fermions and the lower is for bosons. Because the e~^^'^ fac- 
tor is expected to be small for all particles except pions (T ^ 100 MeV) the following 



approximation can be made in such cases 
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and then 
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The 2^(Q° — qj)/2^(Q'') term in Eq.(7) represents classical chemical potential 

Modification of the statistical hadronization model 

The conventional Boltzmann-Gibbs description shown above was modified to allow the 
description of the systems of not-completely-free particles: the correlation "strength" , how- 
ever defined, was introduced with the help of the new non-extensivity parameter and the 
new statistics, which in such case has to be also non- extensive. In the limit of the absence 
of correlations the new description approach the Boltzmann form. 

There could be infinitely many "generalised" statistics which fulfil such requirements. 
Following 'Ockham razor' we should choose the one which is simple and has well defined 
theoretical background. In this present paper we test the possibility, proposed by Tsallis 
[11], based on the modification of the classical entropy definition 

w 
Sbg = -k^ Pi^npi (8) 



by the new one 

S. = kj^^(^l-f:p?j (9) 

with the new parameter q called the non-extensivity parameter. This modification has been 
adopted in other physical applications (see, e.g., [12]). 

Maximisation of the entropy requirement with the total energy constraint 

^^ = Eo (10) 



leads to the probability of realisation of the state i (with energy Ei)gwen by 

P? = ^ [1 - (1 - <l)/T,{E, - Eo)f^'-'^^ , (11) 

where Zg is the normalisation constant related to 2(Q°) of Eq.(l) where the Boltzmann 
terms x could be replaced by the probabilities of the form given in Eq.(ll). 
Equation (11) can be rewritten introducing new symbol, e^, defined as 

(e,)- = [H-(l-g)x]^/(^-^) (12) 

(for completeness, we should mention that in the q=l limit we have, as we should, e^ = e^). 
Then 

p1 - [1 - (1 - Q)/Tq{E, - EoW^^'-"^ ~ e-^'/^' (13) 

and the modified partition function can be written in the form 

2n 2tt 2tt 

states *> n n n 

RESULTS 

We have evaluated Zq functions (and Z = Zi) for a variety of thermodynamical parameter 
values T, V , and strangeness for Q values which cover the production of over 100 hadrons 
of masses below 2 GeV/c^. All decays of short-lived particles were then performed. In the 
Ref. [13] we have shown the effect of exchange the statistical factor from the Boltzmann 
one to the Tsallis modified exponential factor on the p± distribution. We have found that 
the description of transverse momenta with non-extensive statistics works quite well and 
follow the earlier found [14] trends. 

The measurements of some identified particle ratios done mostly by the ALICE Collab- 
oration gives the opportunity to test the statistical model of particle flavours creation in 
the new, higher energy region. The lower energy jet hadronisation results have been anal- 
ysed by the serial of papers by Becattini and others ([15, 16]). It has been shown that the 
micro-canonical Boltzmann description works well for e+e^ from ^/s ?^10 [17] to 91 GeV at 
Ref. [15] and pp and pp interactions up to SPS evergies, y/s ^900 GeV [16]. 

The comparison of the results obtained with Boltzmann (non-extensivity parameter q = 
1) and Tsallis (with the value of g = 1.12 as expected from Ref. [13]) is shown in Fig. 1. We 
can see the clear enhancement of the exponential 'tail'. 
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FIG. 1: Relative particle multiplicities, / = (ni/rich), obtained for the Tsallis statistics (with 
T = 150, V = 100 and q = 1.12) compared with the same ratios for the Boltzmann statistics (with 
T = 150, V = 100 and q = 1). 



The effect of generalisation of the statistics is seen also in the Fig 2. The strong effect 
observed when we compare Boltzmann q parameter 1.0 results with the higher q values is 
clear. The analysis of LHC data on transverse momenta suggests that the appropriate value 
at 7 TeV centre of mass energy is about equal to 1.12. For small changes of q the differences 
are still not very significant. 

We have checked different particle multiplicities dependencies on main thermodynamical 
hadronization parameters T and V. Some results are shown, as an example, in Fig. 3 

The effect of volume V have to be taken with care, because the spacial and temporal 
history of the hadronisation process is not exactly known. It is expected that the canonical 
picture should take into account the multichain idea (e.g., [10, 18]) of decomposition of the 
'hadronic soup' to chain of independently hadronized objects/fireballs. For the Boltzmann 
statistics, by the definition of extensivity, the sum of many hadron sources is equivalent to 
one big source [17]. This is, in general, not the case of Tsallis non-extensive statistics. But 
we can say that the strength of the non-extensivity is still rather small and the effect of 
subdivision of the hadronisation volume does not change the conclusion about the identified 
particles ratios. Another point is related to the effect of the canonical treatment of small 
fireballs which relates to the suppression of strange quark (and diquark, or strange diquark) 
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FIG. 2: Relative particle multiplicities, / = {ni/rich), obtained for Boltzmann and different Tsallis 
generalised statistics (with T = 150, F = 100). 
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FIG. 3: Relative multiplicities, / = (ni/rich), calculated for Tsallis statistics for different temper- 
atures T, (left) and different hadronization volumes V (right). 

what was mentioned already in [1]. We will come back to these additional suppressions later 
on. 

Fig. 3 suggest that ratios of identified particle multiplicities are not sensitives to the 
thermodynamical parameters of the model of fireballs hadronization. It seems that the 
predictions can not be adjusted to the data, at least in the range of possible, reasonable 
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changes, which is quite small. They are as came out from integrations in Eq.(7) and modified 
statistics Eq.(14) with q fixed by p_L distribution fit. The comparison with the data could 
be thus the experimentum crucis of the thermodynamical models in general. Results for 
particle ratios are shown in comparison with data measured by the ALICE Collaboration 
in Tab. I. Thermodynamical parameter values were taken to reproduce roughly the data, 
where it is possible. 

TABLE I: Ratios of particle multiplicities calculated for Boltzmann and Tsallis statistics with spin 
counting states weight factor equal to (2J + 1) and strangeness suppression 7^ = 0.5 acting of 
strange quark contents compared with the measurement results from ALICE. 







Boltzmann 


Tsallis 


Tsallis 


particle ratio 


ALICE measurement 


r=160 
g=1.00 


V=80 
T=170 
g=1.12 


V=80 
T=150 
g=1.15 


p/u 


1.15 ± 0.2 ±0.12^* 


0.985 


0.855 


0.848 


<A/(P + ^) 


0.084 ± 0.013 ± 0.012^* 


0.042 


0.035 


0.033 


K*o/K- 


0.35 ± 0.001 ± 0.04^* 


0.337 


0.466 


0.466 


(j)/K*^ 


0.33 ± 0.004 ± 0.05* 


0.268 


0.215 


0.207 


cP/tt- 


0.014 ± 0.0002 ± 0.002* 


0.0063 


0.0080 0.0077 


<t^/K~ 


0.11 ±0.001 ±0.02* 


0.090 


0.100 


0.097 


w/vrO 


0.6 ±0.1^ 


1.36 


0.861 


0.704 


Q/S 


0.067 ±0.0^^ 


0.068 


0.237 


0.240 


^/<P 


0.04 ± .008" 


0.119 


0.362 


0.403 


ry/vr" 


0.1067 ±0.0259 ±0.0212^ 


0.206 


0.092 


0.081 



'^1 GeV/c < p± < 5 GcV/c [19] 
''full phase space [5] 
''p^>2.5 GeV/c [20] 
'^mi_ -mo> 0.3 GeV [9] 
^pi_ > 1 GeV [5] 
f'pjL > 0.55 GeV/c [8] 



We clearly see problems when comparing the ALICE results with any model predictions, 
whatever parameter set and model is taken. It could be the substantial and general disagree- 
ment (at least in some cases), however, we have the weight factor Wj in Eqs. (5,7) which give 
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FIG. 4: Relative particle multiplicities, / = (ni/rich), obtained for the Tsallis statistics and different 
values of the strangeness suppression parameter 7^ parameter (T = 150, V = 100, q = 1.12). 

us some hope and freedom to fit the model predictions. The expected form of (2 J + 1) is, 
in general, modified since there had been experimentally found the suppression of Ks with 
respect to non-strange particles. The factor 7s was introduced and adjusted for obtain the 
correct value of K/n ratio. The general statement is that the strange phase space is not full 
available for particle production what can be realised by multiplying the partition function 
by given factor for each strange valence quark in the particle in question. The effect is of 
course significant for heavy, strange mesons, but not only. It is shown in Fig. 4 

The strangeness suppression factor is one of the basic parameters in the jet fragmentation 
model introduced by Feynman and developed finally by the Lund group [6]. In the Lund jet 
fragmentation process new hadrons appear as a breaking of the colour field string stretched 
between quarks moving apart by the production of a new pair of quarks (sometimes di- 
quarks). If there is enough energy left, further breaks may occur and eventually only on- 
mass-shell hadrons remains. The creation of new quark-antiquark pair in the Lund model 
is a kind of the quantum tunnelling process, so it is expected that the heavy quark creation 
is suppressed. It is usually assumed that -u : rf : s ~ 1 : 1 : 0.3. 

Additional Wj weight factor is related with the spin stated of newly created hadrons: 
for mesons pseudoscalar and vector states. The eventual suppression here is not defined in 
the Lund fragmentation model. The counting the spin states gives 1 : 3 ratio, but in the 
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JETSET model it is modification according to the 'tunnelling normalisation' and this ratio 
is eventually close to 1 : 1. 

The situation with baryon creation in the Lund model is much more complicated. The 
tunnelling mechanism is also adopted here. We have here the probability of string breakup 
via diquark mode and further combination of quark and diquark. If we take into account 
the pop-corn mechanism of diquark breakups, and lack of general rules, we can have the 
number of parameters to be adjusted to the data comparable with the number of measured 
ratios to be used for the adjustment. The number of parameters describing the production 
of baryons measured with good accuracy in the experiments at LHC is at the moment higher 
then the number of such baryons itself [6] . 

The Lund model and particular JETSET hadronization generator is used also by the 
PHOJET [7] program package for the recent theoretical examination of the LHC data de- 
scription and comparison. Some parameters in PHOJET are different than the default Lund 
model values. 

We first discuss the possibility of introducing the strangeness suppression factor 7^ for 

Wj = (2J + 1) X (7,)^^- . (15) 

The Nj is the 'degree of strangeness' which is in fact not exactly defined. It should be related 
to the contents of the particle j. Three possibilities are rather natural. 

S strangeness of the particle j 

Ug number of strange (or antistrange) valence quarks of j (16) 

Uss number of ss pairs involved to create the particle j 

The difference could be seen in the comparison of K and weights. For direct K they are 
7s, 7s, and 7^, while for direct they are 1, 7^ , and 7^, for the first, the second and the third 
possibility, respectively. The actual situation is more complicated, because of the effects of 
decays heavy resonances, so to see the effect complete calculations have to be performed. 
Some example results are given in the Tab. II. 
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It can be seen that calculated ratios shown in Tab. II that the case {Nj = n^) works well 
for strange mesons and the simple (2 J + 1) factor results not far from measurements. 
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TABLE II: Ratios of particle multiplicities calculated for T = 150, y = 80 and q = 1.12 with 
different strangeness suppression factor definition with the strangeness suppression factor equal to 
75 = 0.5 compared with the measurement results. 



particle ratio 


ALICE measured ratios 


Us 


Uss S 


p/uj 


1.1 ±0.2 ±0.12" 


0.862 


0.870 0.888 


(/./(p + w) 


0.084 ± 0.013 ± 0.012^* 


0.0360 0.0718 0.142 


K*0/K' 


0.35 ± 0.001 ± 0.04^* 


0.462 


0.423 0.382 


ct>/K*^ 


0.33 ± 0.004 ± 0.05^ 


0.219 


0.440 0.751 


0/vr- 


0.014 ±0.0002 ±0.002^ 


0.0080 0.0160 0.0301 


(j^/K- 


0.11 ±0.001 ±0.02^ 


0.101 


0.186 0.287 


uj/tt^ 


0.6±0.P 


0.894 


0.893 0.870 


0/H 


0.067 ±0.01^^ 


0.233 


0.219 0.233 


n/<i) 


0.04 ± .008" 


0.329 


0.135 0.082 


ry/vrO 


0.1067 ±0.0259 ±0.0212-^ 


0.100 


0.102 0.103 



Further 'tuning' can involve the value of the strangeness suppression factor 7^. we have 
checked tree values which have been used in the literature: 0.5 originally proposed by Feyn- 
man in the jet fragmentation model [2] and still of use [10] , 2/3 used successfully by Becattini 
[21] for the e^e~ data and pp results from SPS, and 7^ = 1 as the limit of no strangeness 
suppression. The results are shown in Tab. Ill 



The general agreement, at least for mesons, is seen for 7^ k. 0.5 — 0.66. 

The discrepancy exists still in ratios involved baryons. As it has been said there is a great 
degree of freedom to modify baryon multiplicities. The diquark suppression factor 7^^ is the 
one possibility which we used and another factor was introduced specially for Vt baryons 7^5 
related to creation of double strange diquark. 

In the Tab. IV we present results for the 7^^ and 7^5 equal to 0.5. 



Because of relatively limited amount of data we do not wish at the moment to go further 
with 'tuning' the parameters (7^, 7gg and 7^^). 
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TABLE III: Ratios of particle multiplicities calculated for T = 150, y = 80 and q = 1.12 with 
different strangeness suppression (for Us) factor values and spin counting states exactly equal to 
(2 J + 1) compared with the measurement results. 



particle ratio 


ALICE results 


7s = 1 


7, = 0.66 


7s = 0.5 


p/u 


1.15 ±0.2 ±0.12'' 


0.882 


0.867 


0.861 


0/(p + ^) 


0.084 ± 0.013 ± 0.012'=' 


0.137 


0.0618 


0.0360 


K*0/K- 


0.35 ± 0.001 ± 0.04^ 


0.429 


0.453 


0.462 


HK*^ 


0.33 ±0.004 ±0.05^ 


0.433 


0.289 


0.219 


4>h~ 


0.014 ± 0.0002 ± 0.002^ 


0.0234 


0.0126 


0.00802 


4>IK- 


0.11 ±0.001 ±0.02^ 


0.186 


0.131 


0.101 


uj/tt'^ 


0.6±0.r 


0.787 


0.857 


0.894 


n/E 


0.067 ± O.Ol'^ 


0.443 


0.303 


0.233 


n/cp 


0.04 ± .008^ 


0.585 


0.416 


0.329 


7]/tT^ 


0.1067 ±0.0259 ±0.0212^ 


0.084 


0.094 


0.161 



The ALICE Collaboration published data showing the particle ratio as a function of 
transverse momentum of the particles. Taking into account that the modification of the 
statistics of the multiparticle production process was developed primary for the transverse 
momentum description this kind of data could be valuable to verify the model. Comparison 
of our final model prediction and the data are shown in Fig. 5. 

The Boltzmann statistics results of the fit with parameters of Ref. [16] are also given for 
comparison. As it is seen, the standard statistics does not work very well for the LHC ALICE 
data shown in Fig. 5, while the modified, Tsallis statistics with our chosen suppression factors 
reproduces the data quite better. 

CONCLUSIONS 

We have shown that the introduction of non-extensive statistics to the thermodynamical 
theory of multiparticle production in hadronic collisions can give a consistent description 
of identified particle multiplicity ratios. The thermodynamical model parameters found 
analysing the transverse momentum distributions measured at 7 TeV could be used without 



13 



TABLE IV: Ratios of particle multiplicities calculated for T = 150, y = 80 and q = 1.12 with 
different strangeness suppression factor 0.5 and 0.66 (and jqq = ^ss = 0.5) compared with the 
measurement results. 



particle ratio 


ALICE results 


7. = 0.5 


-is = 0.66 


p/uj 


1.15 ± 0.2 ±0.12^* 


0.862 


0.867 


cl)/{p + uj) 


0.084 ± 0.013 ± 0.012°^ 


0.036 


0.0620 


K*o/K- 


0.35 ± 0.001 ± 0.04^* 


0462 


0.453 


ct>/K*^ 


0.33 ± 0.004 ± 0.05^ 


0.226 


0.302 


0/vr- 


0.014 ± 0.0002 ± 0.002^ 


0.0091 


0.0147 


<t^lK- 


0.11 ± 0.001 ± 0.02^* 


0.104 


0.136 


L^/vrO 


0.6 ±0.1^ 


0.907 


0.872 


0/H 


0.067 iO.Ol'^ 


0.0797 


0.103 


n/cj) 


0.04 ± .008^ 


0.0363 


0.0447 


r^/vr" 


0.1067 ±0.0259 ±0.0212^ 


0.111 


0.105 



readjustment for identified particle multiplicity data description. The standard strangeness 
suppression factor 7 of about 0,5-0.6 with additional suppression of diquark and strange 
diquark phase space works well for average as well as p_[_ dependence of particle ratios. 

Proposed treatment of hadronization can put some new light on the particle creation 
physics and it is somehow complementary to the jet type non-perturbative QCD models. 
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FIG. 5: Identified particles ratios as a function of transverse momentum (transverse mass for 0/H) 
for T = 150, Y = 80, q = 1.12 and strangeness suppression 7^ = 0, 5 acting of strange quarks (or 
antiquarks) contents of the particle in comparison with ALICE data. 
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